1 
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o : 

We study qqqq-type four-quark (4Q) systems in SU(3)c anisotropic quenched lattice 

1— ( I QCD, using the 0(a)-improved Wilson (clover) fermion at /3 = 5.75 on 12^ x 96 with renor- 

^ . malized anisotropy as /at = 4. For comparison, we first investigate the lowest qq scalar 

' ' ' meson from the connected diagram and find its large mass of about 1.32GeV after chiral 

CO ' extrapolation, and thus the lowest qq scalar meson corresponds to /o(1370). We investigate 

^SJ , the lowest 4Q state in the spatially periodic boundary condition, and find that it is just a 

■ two-pion scattering state, as is expected. To examine spatially-localized 4Q resonances, we 
' ^ / ' use the Hybrid Boundary Condition (HBC) method, where anti-periodic and periodic bound- 
ed ' ary conditions are imposed on quarks and antiquarks, respectively. By applying HBC on a 

, finite-volume lattice, the threshold of the two-meson scattering state is raised up, while the 

^ • mass of a compact 4Q resonance is almost unchanged. In HBC, the lowest 4Q state appears 

^ ' slightly below the two-meson threshold. To clarify the nature of the 4Q system, we apply the 

I Maximum Entropy Method (MEM) for the 4Q correlator and obtain the spectral function 

I ' , of the 4Q system. From the combination analysis of MEM with HBC, we finally conclude 
that the 4Q system appears as a two-pion scattering state and there is no spatially-localized 

' 4Q resonance in the quark- mass region of nis < niq < 2ms. 

a^ : 

O ■ 81. Introduction 

cn ■ ^ 
<^ ■ 

; There are five 0++ isoscalar mesons below 2GeV: /o(400-1200), /o(980), /o(1370), 

Q ■ /o(1500) and /o(1710). Among them, /o(1370) is assigned to be the lowest qq scalar 

. meson in the quark model^ since the lowest qq scalar meson is ^Pq in the quark 

^p. \ model and it turns to be rather heavy. /o(1500) and /o(1710) are expected to be 

■ the lowest scalar gluebalP''^ or an ss scalar meson.'S Then, what are the two light 
^ , scalar mesons, /o(400-1200) and /o(980)? This is the "scalar meson puzzle", which is 

^ \ unsolved even at present. As a possible answer, in 1977, JafFe proposed tetra-quark 

(qqqq) assignmeniP for low- lying scalar mesons such as /o(980) and ao(980). 

The tetra-quark is also interesting in terms of the recent experimental discoveries 
of charmed tetra-quark candidates, e.g., Di(2317), D+i(2460), X(3872) and Y(4260), 
which are expected to be tetra-quark statea^''^ from the consideration of their mass, 
narrow decay width and decay pattern. In fact, simple quark-model assignment 
for D^(2317) is cs{^Pq), but the P-wave cs state is much heavier than 2.32GeV.'3 
Furthermore, there is a puzzle in the decay mode of D^(2317): CLEO experiments 
show that the dominant decay mode is Dj^vr*^ and its contribution is much larger 

than the radiative decay as ^Ij^+^j^^ p| < 0.059. If D^(2317) is isoscalar such as cs, 

this decay pattern implies that the main decay mode is isospin breaking process of 
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O(e^) in QED, which is rather anomalous. To explain this decay pattern, Terasaki 
proposed 1 = 1 possibility of D^(2317)|^ which means the tetra-quark picture for 
D^(2317) and predicts a doubly-charged meson D^^"^ as the isospin partner. 

In this way, multi-quark systems {qqqq, qqqqq, • • • ) have been studied as an 
interesting subject, and, recently, they are also investigated in lattice QCD^''^''^''^ 
as well as ordinary qq and qqq systems .'^''^''^ In this paper, we investigate the 
four-quark (4Q) system and examine the scalar tetra-quark resonance as fo{udud) by 
treating only the connected diagram in quenched anisotropic lattice QCD. Note that 
the tetra-quark D^"*" {cusd) corresponds to the connected diagram of /o (uudd) in the 
unrealistic idealized SU(4)j limit of mu,d=ms='mc. Then, the lattice QCD study of 
the connected diagram of fo{uudd) means that of the manifestly exotic tetra-quark 
D++(cnsJ) with C = +1, S = +1, I3 = +1 in the idealized SU(4)/ limit.l^ 

§2. Anisotropic Lattice QCD and Hybrid Boundary Condition Method 

To examine exotic 4Q states as fo{udud), we adopt a diquark-type interpolating 
field, O = Eabc {uaClbdh) X Edec {udC'y^d^) , which is expected to have a small overlap 
with two-meson scattering states. In lattice QCD, we calculate and analyze the 
temporal 4Q correlator G{t) = y x)O^(0, 0)), where the total momentum 

of the 4Q system is projected to be zero. We impose the Dirichlet boundary condition 
for quarks in the temporal direction to avoid the backward meson propagation.^^' 

To get detailed information on the temporal behavior of the 4Q correlator G{t), 
we adopt anisotropic lattice QCD with the plaquette action and the 0(a)-improved 
Wilson fermion at /? = 2Nc/g^ = 5.75 on 12^ x 96 with renormalized anisotropy 
tts/at = 4. We use 500-1000 gauge configurations. The temporal lattice spacing is 
at — 0.045fm, and the spatial lattice size is L = 12as — 2.15fm. We take fi;=0.1210, 
0.1220, 0.1230 and 0.1240, which cover the quark mass region of < ruq < 2ms. 

To examine low- lying 4Q resonances, we use the "Hybrid Boundary Condition" 
(HBC),I3''S'II3 where we impose the anti-periodic boundary condition (APBC) for 
quarks (it,(i) and the periodic boundary condition (PBC) for antiquarks (n,d), as 
shown in Table I. By applying HBC on a finite lattice with L^, the two- meson thresh- 
old is raised up, while the mass of a compact 4Q resonance is almost unchanged. 
Then, the 4Q resonance may become visible as a low-lying state in HBC, if it exists. 

For a compact qqqq resonance, since it contains even number of quarks, it obeys 
PBC even in HBC, and therefore its energy in HBC is almost the same as that in 
PBC.^' For a two-meson scattering state, both mesons obey APBC and have non- 
zero relative momentum pmin = ±X' ^x)' ^'^'j IPminl = \/3vr/L in HBC. In fact, 
the two-meson threshold m\ + 777.2 ("^ij "^2'- meson masses) in PBC is raised up in 

HBC as ^th - Efc=i,2 +fmm with |pmin| = V^TT / L ~ 0.5GeV for L ~ 2.15fm. 



Table I. Hybrid Boundary Condition (HBC) to raise the threshold of two-meson scattering states. 
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§3. EfTective-Mass Analysis for Four-Quark Systems 

We perform the standard effective-mass analysis for the low-lying 4Q state."^ 
For comparison, we first calculate the ordinary qq scalar-meson mass from the 
connected diagram in anisotropic quenched lattice QCD. We show in Fig. 1(a) the 
lowest qq scalar meson mass plotted against m^. Our quenched lattice QCD indi- 
cates that the lightest qq scalar meson has a large mass as Mqq ~ 1.310(55)GeV 
after the chiral extrapolation. Thus, the lightest qq scalar meson is identified with 
the multiplet of /o(1370) and ao(1450),'^ which is consistent with the quark- model 
assignment for scalar mesons. (If disconnected diagrams are dropped off, isoscalar 
and isovector qq mesons degenerate in quenched QCD, e.g., nip = m^, mf^ = ma^.) 




m/ [GeV^] [GeV^] m/ [GeV^] 



Fig. 1. (a) The lowest mass Mqq of the qq scalar meson plotted against m^. The qq scalar meson 
mass Mqq is extracted from the connected diagram in quenched lattice QCD. (b) The lowest 
mass M4Q of the scalar 4Q (udud) system against in PBC. The curve denotes the two-pion 
threshold, (c) The lowest mass M4Q of the scalar 4Q system against in HBC. Open (closed) 
symbols are the direct lattice QCD data (chiral extrapolated values). The curve denotes the 
raised two-pion threshold Eth in HBC. Triangles denote the qq scalar-meson mass for comparison. 

We calculate the lowest energy of the scalar 4Q (udud) system in PBC and in 
HBC, as shown in Figs.l (b) and (c). In PBC, the lattice QCD data coincide with the 
two-pion threshold curve, and hence the lowest 4Q state is just a two-pion scattering 
state. In HBC, the 4Q state appears slightly below the raised two-pion threshold 
i?th by 40-75MeV and seems to behave linearly in rather than the two-pion 
threshold. Then, there are two possible interpretations for the lowest 4Q state in 
HBC: i) If the slight difference from the two-pion threshold can be explained by the 
two-pion interaction and the finite-box effect in HBC, the lowest 4Q state in HBC is 
also regarded as a two-pion scattering state; ii) From the slight difference from the 
two-pion threshold and different chiral behavior, one may expect that the 4Q state 
in HBC differs from the two-pion scattering state. If a linear chiral extrapolation is 
adopted, the 4Q mass is estimated as m^Q ~ 1.174(18)GeV in the chiral limit.'^ 

§4. MEM Analysis for Four-Quark Systems 

To clarify whether the 4Q resonance exists or not, we perform the Maximum 
Entropy Method (MEM) analysis for the scalar 4Q system in PBC and in HBC, and 
obtain the spectral function A{uj) of the scalar 4Q state from the 4Q correlator G{t) 
with K=0.1240, as shown in Fig. 2. Here, we calculate the default function m{uj) of the 
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Fig. 2. The dimensionless spectral function p{uj) = A{u})/uj^ obtained from the point-source point- 
sink 4Q correlator G(t) for k=0.1240 (i.e., m, ~ m^) in PBC (left) and HBC (right). 

4Q system with the lowest-order perturbative QCD: m(a;) = 2^r{t)r{6)TT'^ ^^ ~ "^o"^^) 
which gives the asymptotic form of A(uj) in the high-energy region. 

In both cases of PBC and HBC, the lowest peak around l.SGeV exhibits a large 
strength. Note however that, in the finite-size lattice, the momentum of the pion 
is discretized, and all the scattering states are observed as resonance- like states in 
lattice QCD. By changing the spatial boundary condition from PBC to HBC, all the 
low-lying peaks in the spectral function are significantly shifted to the high-energy 
side. Since each low-lying peak is strongly affected by the spatial boundary condition, 
it is not a spatially-localized resonance but a spatially-extended two-pion scattering 
state, i.e., a discretized continuum state. (If there is a low- lying pole whose position 
is unchanged between PBC and HBC, it is to be a spatially- localized resonance.'^) 

Thus, from the combination analysis of MEM with HBC, we finally conclude 
that the 4Q system appears as a two-pion scattering state and there is no spatially- 
localized 4Q resonance in the quark-mass region of < niq < 2ms. 
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